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Sigma pole position and errors of a once and twice subtracted dispersive 
analysis of pi-pi scattering data 
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We show how the new precise data on kaon decays together with forward dispersion relations, sum rules and 
once- and twice-subtracted Roy's equations allow for a precise determination of the sigma meson pole position. 
We present a comparison and a study of the different sources of uncertainties when using either once- or twice- 
subtracted Roy's equations to analyze the data. Finally we present a preliminary determination of the a pole 
from the constrained dispersive data analysis. 



1. Introduction 

Roy's equations (RE), based on twice- 
subtracted dispersion relations and crossing sym- 
metry conditions for tttt — > tttt amplitudes were 
obtained in 1971 [1]. In recent years, these equa- 
tions have been used either to obtain predictions 
for low energy tttt scattering, sometimes using 
Chiral Perturbation Theory (ChPT) or to 

test ChPT [41516] , as well as to solve old data am- 
biguities [7]. Roy eqs. are relevant for the sigma 
pole, whose position has also been predicted very 
precisely with the help of ChPT [8]. 

Our group |5|6j has also used Roy eqs. with 
Forward Dispersion Relations (FDR) to obtain 
a precise determination of tttt scattering ampli- 
tudes from data consistent with analyticity, uni- 
tarity and crossing. On purpose, we have not 
included ChPT constraints, so that we can use 
our results as tests of the ChPT predictions. Un- 
fortunately, the large experimental error of the 
scattering length Oq of the isospin 2 scalar partial 
wave, becomes a very large error for the sigma 
pole determination using RE. For this reason, 
a new set of once-subtracted RE, called GKPY 
eqs. for brevity, have been derived. Both the RE 
and GKPY equations provide analytic extensions 



for the calculation of poles in the complex plane. 
Here we present preliminary results for the a pole 
in the SO wave obtained from GKPY eqs. and 
analyse the different sources of errors comparing 
their size with those from RE. 

2. Roy versus GKPY Equations 

In practice Roy's equations relate real parts 
(output) of the 5*0, P and 5*2 partial wave ampli- 
tudes //(s) below some energy scale Smax, with 
the imaginary parts of all other partial waves 
and high energy parametrizations (input). Both 
GKPY eqs. and RE can be split into three terms: 

i?e//(s) = STiis) + KTl{s) + DT,'{s). 

The "subtraction terms" ST{s) are polynomi- 
als whose coefficients are linear combinations of 
the SO and S2 scattering lengths. For standard, 
twice subtracted, Roy's eqs., ST are first degree 
polynomials, thus growing quadratically with en- 
ergy, whereas for GKPY they are just constants. 
Hence, the uncertainty due to the poor experi- 
mental knowledge of the isospin 2 scalar scatter- 
ing length, Uq becomes a very large source of error 
at high energies for standard Roy's eqs. 

The "kernel terms" KT{s) parametrize the 
contribution of the 50, P and S2 waves {£ — 0, 1) 
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up to a given energy s^ax, whereas the "driving 
terms" DT{s) refer to all waves above s^^ax a-nd 
also to £ > 1 partial waves from threshold. In our 
case Y^Smaa: = 1420 MeV and for higher energies 
we use Regge parametrizations. Both the kernel 
and driving terms are integrals of the kind: 

/ ds'i^//'(s,.')Im//;(5'). (1) 

i'=o e •' 

The integration kernels KH (s, s') are different 
for RE and GKPY eqs. The important fact being 
that, for high s', they decrease as ~ 1/ s'^ for RE, 
but as ~ l/s'2 for GKPY eqs. Hence GKPY eqs. 
are more sensitive to the high energy input. 

By comparing Fig. [1] for RE with Fig. [2] for 
GKPY eqs. very remarkable differences (espe- 
cially for the 50 wave) can be observed in the 
relative sizes of the ST and KT terms. Please 
note the different scales. The curves represent 
the real part (output) of the amplitudes, i.e., 
gi/2 ^ gjjj^ 5/2fc, where 5 and rj are the phase shifts 
and elasticities for the preliminary amplitudes of 
the Constrained Data Fit described in the talk by 
J.R. Pelaez in this conference. This fit describes 
data and has been constrained to satisfy FDR, 
RE. GKPY eqs. and some crossing sum rules. 

Note that, for RE, the ST and KT terms, 
which are huge, suffer a strong cancellation. Ac- 
tually, for sufficiently large energy, both the ST 
and KT terms are much larger than the unitarity 
bound \Ret\ < s^/'^/2k ^ 1, which is only satis- 
fied by the real part of the total amplitude after 
their strong cancellation. In contrast, GKPY eqs. 
are dominated by KT terms, which are always 
smaller than one. Remarkably, the DT are still 
relatively small compared with the KT despite 
there is one less subtraction in GKPY eqs. This 
means that the effect of high energy input, whose 
knowledge is less detailed, is well under control. 

The error bands depicted in Figs. [Hand [2] are 
generated using a Monte Carlo Gaussian sam- 
pling of all parameters of the Constrained Fit 
to Data, (varied within 3 standard deviations). 
Uncertainties are just the widths of independent 
Gaussian fits to the left and right sides of the 
ST, KT and DT distributions. As seen in Fig. [TJ 
for sufficiently high energies, the linear s depen- 



dence of the ST terms makes their uncertainty 
large and dominant when using RE for the SO and 
52 waves. This does not happen for GKPY eqs. 
whose errors, given the same input, are smaller, 
above s ~ 400 MeV, than those from RE (see 
also Fig. 1 in J.R. Pelaez's talk in this confer- 
ence). For this reason we are now implementing 
the GKPY eqs. together with FDR and standard 
Roy. eqs. not only to constraint the amplitude 
above ~ 400 MeV, but also to extend analytically 
the amplitudes to determine the a pole position. 

3. Position of the a pole 

The mass and width of the a or /o(600) meson 
quoted in the Particle Data Table are very widely 
spread [9] 

AU - « Y '''^^^ ~ ^'^'^^^ ~ " 500)MeV.(2) 

The main reason of these uncertainties is that tttt 
scattering data are few and sometimes contradic- 
tory. Moreover, all the quoted theoretical models 
are not equally reliable, and even less so when 
extending the amplitude to the complex plane. 
Thus the position of the sigma pole in various 
models differ significantly [S], although, with a 
couple of exceptions, they tend to agree around 
M„ - « (400 - 500) - i(220 - 300) MeV. 

The recent data from E865 collaboration at 
Brookhaven [TOj and from NA48/2 [HJ provide 
us with new and very precise information on the 
TTTT scattering at low energies. Thanks to these 
new data we are able to construct, with our Con- 
strained Fits to Data, a very reliable description 
for the 50 wave especially near the tttt threshold 
(see J.R. Pelaez talk in this conference). 

With those precise data parametrizations, we 
can now use either the Roy eqs. or the GKPY 
eqs. to extend the partial waves analytically to 
the complex plane and look for poles in the second 
sheet of the S-matrix. As it well known, a pole 
on the second Riemann sheet (unphysical sheet) 
is associated with a zero on the first — the physical 
one. Therefore, as usual, we look for zeroes of the 
physical sheet of the 5-matrix, 

50(s) = l + 2,Vl-4m2/s/0(s). (3) 
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Depending on whether we use Roy or GKPY eqs. 
we find a different accuracy in our results, namely 

^ = - i 257tll ^eV (RE) (4) 

^ = AGltlti - i 255 ± 16 MeV (GKPY) (5) 

from the GKPY equations. These values are in 
good agreement with each other. Note that the 
errors have been calculated in the same way as 
for the output amplitudes in Sect. 1, i.e. from 
a Monte Carlo Gaussian sampling. Except for 
those of the mass from RE, which are roughly 
25% smaller, they are almost identical to the non- 
Gaussian symmetric errors that we have been us- 
ing in previous works |12j 

(461±14)-i (255±16)MeV (GKPY)(6) 

It is reassuring that our error estimates remain 
stable calculating them in these two rather differ- 
ent ways. 

This preliminary determination is quite consis- 
tent with other recent and precise results in the 
literature: On the one hand, both the mass and 
width lie less than 1.25 standard deviations from 
the prediction of twice-subtracted Roy's equa- 
tions combined with ChPT results for the scatter- 
ing lengths [8]: = 4411^*^ - i 272l?4 gMeV. 
On the other hand our preliminary mass deter- 
mination above is also slightly beyond one stan- 
dard deviation from the pole in our simple fit 
|13j of a conformal expansion to low energy data 
^ = (484 ±l7)-i (255 ± 10)MeV, but agrees 
remarkably well on the width. 

4. Conclusions 

The GKPY equations - Roy-like dispersion re- 
lations with one subtraction for the tttt ampli- 
tudes - provide stringent constraints for disper- 
sive analysis of experimental data. The main ad- 
vantage of GKPY eqs. is that, for the same input, 
in the 0.4 GeV < y/s < 1.1 GeV region they have 
significantly smaller errors than standard Roy. 
eqs. Hence, they provide better accuracy tests 
and analytic extensions of the amplitudes in that 
region. In particular, using just a data analysis 
consistent within errors with Forward Dispersion 
Relations, Roy eqs. and GKPY eqs. (and no 



ChPT input), we have presented here the follow- 
ing preliminary but very precise determination of 
the cr pole position: 

= 46lii5 j - i 255 ± 16 MeV (7) 

The details of GKPY eqs. together with a full 
constrained dispersive analysis the data, extend- 
ing both GKPY and Roy eqs. up to 1.1 GeV and 
some slight improvements on the partial waves is 
about to be finished. A final number for the a 
pole should follow relatively easy once that data 
analysis is completed. 
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Figure 1. Subtraction (ST), kernel {KT) and 
driving {DT) terms for the 50, P and S2 
waves from the twice-subtracted Roy's equations. 
Dashed bands denote the errors of these terms. 
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Figure 2. As in Fig. [1] but for GKPY equations. 



